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LONG-TERM  GOALS 

We  are  part  of  a  multi-institutional  research  team  funded  by  the  ONR-sponsored  Radiance  in  a 
Dynamic  Ocean  (RaDyO)  program.  The  primary  research  goals  of  the  program  are  to  (1)  examine 
time-dependent  oceanic  radiance  distribution  in  relation  to  dynamic  surface  boundary  layer  (SBL) 
processes;  (2)  construct  a  radiance-based  SBL  model;  (3)  validate  the  model  with  field  observations; 
and  (4)  investigate  the  feasibility  of  inverting  the  model  to  yield  SBL  conditions.  Our  goals  are  to 
contribute  innovative  measurements,  analyses  and  models  of  the  sea  surface  roughness  at  length  scales 
as  small  as  a  millimeter.  This  characterization  includes  microscale  and  whitecap  breaking  waves. 

The  members  of  the  research  team  are 

Michael  Banner,  School  of  Mathematics,  UNSW,  Sydney,  Australia 

Johannes  Gemmrich,  Physics  and  Astronomy,  UVic,  Victoria,  Canada 

Russel  Morison,  School  of  Mathematics,  UNSW,  Sydney,  Australia 

Howard  Schultz,  Computer  Vision  Laboratory,  Computer  Science  Dept,  U.  Mass.,  Mass 

Christopher  Zappa,  Lamont  Doherty  Earth  Observatory,  Palisades,  NY 

OBJECTIVES 

Nonlinear  interfacial  roughness  elements  -  sharp  crested  waves,  breaking  waves  as  well  as  the  foam, 
subsurface  bubbles  and  spray  they  produce,  contribute  substantially  to  the  distortion  of  the  optical 
transmission  through  the  air-sea  interface.  These  common  surface  roughness  features  occur  on  a  wide 
range  of  length  scales,  from  the  dominant  sea  state  down  to  capillary  waves.  Wave  breaking  signatures 
range  from  large  whitecaps  with  their  residual  passive  foam,  down  to  the  ubiquitous  centimeter  scale 
microscale  breakers  that  do  not  entrain  air.  There  is  substantial  complexity  in  the  local  wind-driven  sea 
surface  roughness  microstructure,  including  very  steep  nonlinear  wavelets  and  breakers.  Traditional 
descriptors  of  sea  surface  roughness  are  scale-integrated  statistical  properties,  such  as  significant  wave 
height,  mean  squared  slope  [e.g.,  Cox  and Munk,  1954]  and  breaking  probability  [e.g.,  Holthuijsen  and 
Herbers,  1986].  Subsequently,  spectral  characterisations  of  wave  height,  slope  and  curvature  have  been 
measured,  providing  a  scale  resolution  into  Fourier  modes  for  these  geometrical  sea  roughness 
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parameters.  More  recently,  measurements  of  whitecap  crest  length  spectral  density  [e.g.,  Gemmrich  et 
al.,  2008;  Phillips  et  al.,  2001]  and  microscale  breaker  crest  length  spectral  density  [e.g.,  Jessup  and 
Phadnis,  2005]  have  been  reported. 

Our  effort  seeks  to  provide  a  more  comprehensive  description  of  the  physical  and  optical  roughness  of 
the  sea  surface.  We  will  achieve  this  through  the  analysis  of  our  suite  of  comprehensive  sea  surface 
roughness  observational  measurements  within  the  RADYO  field  program.  These  measurements  are 
designed  to  provide  optimal  coverage  of  fundamental  optical  distortion  processes  associated  with  the 
air-sea  interface.  In  our  data  analysis,  and  complementary  collaborative  effort  with  RaDyO  modelers, 
we  are  investigating  both  spectral  and  phase-resolved  perspectives.  These  will  allow  refining  the 
representation  of  surface  wave  distortion  in  present  air-sea  interfacial  optical  transmission  models. 

APPROACH 

We  build  substantially  on  our  accumulated  expertise  in  sea  surface  processes  and  air-sea  interaction. 
We  are  working  within  the  larger  team  (listed  above)  measuring  and  characterizing  the  surface 
roughness.  This  team  is  contributing  the  following  components  to  the  primary  sea  surface  roughness 
data  gathering  effort  in  RaDyO: 

•  polarization  camera  measurements  of  the  sea  surface  slope  topography,  down  to  capillary  wave 
scales,  of  an  approximately  lmx  lm  patch  of  the  sea  surface  (see  Figure  1),  captured  at  video 
rates.  [Schultz,  Zappa] 

•  co-located  and  synchronous  orthogonal  75  Hz  linear  scanning  laser  altimeter  data  to  provide 
spatio-temporal  properties  of  the  wave  height  field  (resolved  to  0(0. 5m)  wavelengths)  [Banner, 
Morison] 

•  high  resolution  video  imagery  to  record  whitecap  data  from  two  cameras,  close  range  and  broad 
field  [Gemmrich] 

•  fast  response,  infrared  imagery  to  quantify  properties  of  the  microscale  breakers,  and  surface  layer 
kinematics  and  vorticity  [Zappa] 

•  air-sea  flux  package  including  sonic  anemometer  to  characterize  the  near-surface  wind  speed  and 
wind  stress  [Zappa] 


The  team’s  envisaged  data  analysis  effort  includes:  detailed  analyses  of  the  slope  field  topography, 
including  mean  square  slope,  skewness  and  kurtosis;  laser  altimeter  wave  height  and  large  scale  wave 
slope  data;  statistical  distributions  of  whitecap  crest  length  density  in  different  scale  bands  of 
propagation  speed  and  similarly  for  the  microscale  breakers,  as  functions  of  the  wind  speed/stress  and 
the  underlying  dominant  sea  state.  Our  contributions  to  the  modeling  effort  will  focus  on  using  RaDyO 
data  to  refine  the  sea  surface  roughness  transfer  function.  This  includes  the  representation  of 
nonlinearity  and  breaking  surface  wave  effects  including  bubbles,  passive  foam,  active  whitecap  cover 
and  spray,  as  well  as  micro-breakers. 

WORK  COMPLETED 

Our  effort  in  FY13  comprised  further  collaborative  analysis  of  the  sea  surface  roughness 
measurements  gathered  from  FLIP  during  the  RaDyO  field  experiments  in  the  Santa  Barbara  channel 
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during  Sept.  5-27,  2008  and  Hawaii  during  Aug.  23-Sept.vl6,  2009.  We  participated  actively  in 
extending  our  analysis  and  validation  of  the  polarimetric  and  wave  breaking  data.  We  submitted  two 
publications.  One  is  in  press  [Gemmrich  et  al.,  2013]  and  one  is  in  revision  [Banner  et  a/.,  2013]. 

RESULTS 

To  provide  the  context  of  our  new  results  for  FY13,  during  the  SBC  experiment,  diurnal  processes 
were  the  origin  of  most  of  the  variability  in  the  atmospheric  and  oceanic  forcing.  The  momentum, 
sensible  heat  and  latent  heat  fluxes  all  show  a  diurnal  cycle  similar  to  the  wind  speed.  In  strong 
contrast,  in  the  HI  experiment,  processes  were  driven  by  the  persistent  easterly  trade  winds.  This 
overview  of  the  background  air-sea  flux  conditions  was  obtained  from  time  series  of  momentum  (wind 
stress),  sensible  heat,  and  water  vapor  direct  covariance  fluxes  during  the  Santa  Barbara  Channel 
(hereafter  SBC)  and  central  Pacific  Ocean  studies  of  Hawaii  (hereafter  HI).  Full  details  are  provided  in 
Zappa  et  al.  [2012]  and  Dickey  et  al.  [2012]. 

In  response  to  the  growing  need  for  robust  validation  data  for  the  Phillips  [1985]  breaking  wave 
spectral  framework,  we  contribute  new  field  results  observed  from  R/P  FLIP  for  the  breaking  crest 
length  distributions  during  two  different  wind-wave  conditions,  and  breaking  strength  during  one 
wind-wave  condition.  According  to  Phillips  [1985],  the  scale  of  each  breaking  front  wave  is  associated 
with  the  phase  speed  of  the  underlying  wave  that  is  breaking,  and  the  associated  width  of  the  breaking 
front  contributes  at  that  scale  to  the  distribution  of  mean  breaking  crest  width/unit  sea  surface  area, 
A(c).  Here  we  characterize  each  breaking  front  velocity  by  its  initial  observed  velocity  Cb,  where  Cb=ae, 
where  a~0.8,  and  define  the  azimuthally  integrated  crest  length  distribution  of  measured  A(Cb).  Then, 
according  to  Phillips  [1985]  breaking  wave  spectral  framework,  the  spectral  distribution  of  wave 
energy  dissipation  rate  e(Cb)  is  given  by  e(Cb)  =  b(Cb)  p  g’1  Cb’5  A(cb),  where  b(Cb)  is  the  scale 
dependent  breaking  strength  coefficient  that  links  the  breaker  front  kinematics  to  the  energetics  and 
dynamics  of  the  breaking  waves.  The  water  density  is  p  and  g  is  the  gravitational  acceleration. 

Breaking  crest  length  distributions 

The  first  experiment  in  the  SBC  had  developing  seas  and  the  second  experiment  near  HI  had  mature 
seas.  These  are  amongst  the  first  experiments  to  use  dissipation  rate  measurements  probing  up  into  the 
breaking  crest  together  with  simultaneous  measurements  of  breaking  crest  length  distributions.  Of 
particular  interest  is  that  the  dynamic  range  of  the  wind  stress  is  larger  during  the  SBC  experiment  than 
the  Hawaii  experiment.  Together  these  two  experiments  provide  wave  breaking  measurements  over  an 
interesting  dynamic  range  of  wind  speeds  spanning  light  and  variable  to  strong. 

Results  on  the  breaking  crest  length  distribution  A(Cb)  for  the  two  experiments  are  given  in  Fig.  2.  A 
key  feature  of  all  observed  A(Cb)  distributions  at  both  experimental  sites  is  a  maximum  at  short-to- 
intermediate  wave  scales.  The  location  of  this  maximum,  which  specifies  the  scale  of  the  largest 
contribution  of  breaking  crests,  increases  with  wave  age,  however  there  is  no  clear  dependence  on 
wave  age  of  the  magnitude  of  this  maximum.  Wind  speed  and  wave  field  parameters  covered  a  much 
broader  range  during  the  Santa  Barbara  Channel  experiment  compared  to  the  experiment  in  the  central 
Pacific  Ocean  south  of  Hawaii.  These  differences  in  the  dynamic  range  of  forcing  parameters  are 
directly  reflected  in  the  dynamic  range  of  the  breaking  crest  length  distributions.  At  intermediate 
scales,  i.e.  3  m/s  <  Cb  <  5  m/s,  the  A(cb)  distributions  span  almost  two  orders  of  magnitude  in  SBC,  but 
only  a  factor  0(5)  in  HI.  The  maximum  levels  of  A(cb)  are  comparable  between  the  two  experiments, 
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but  the  overall  mean  level  of  A(cb)  is  larger  in  HI.  It  is  seen  that  the  slope  of  the  A(cb)  curves  for  higher 
Cb  values  fluctuates  about  the  Phillips  [1985]  canonical  fonn  A(Cb)  ~  Cb’6,  which  is  discussed  further 
below. 


Figure  2.  Distributions  of  spectral  breaking  crest  length  per  unit  area,  A(cb),  against  observed 
breaker  speed  c*  during  the  experiments  in  Santa  Barbara  Channel  (a-c)  and  Hawaii  (d-f).  Colors 
represent  wave  age  cfu*.  The  dashed  line  indicates  the  m  =  -6  slope  predicted  by  Phillips  [1985]. 
For  each  experiment,  data  runs  are  split  into  3  equal-sized  sub-panels,  covering  the  wave  age  sub¬ 
ranges  during  the  entire  data  sets. 


There  are  some  noteworthy  differences  in  the  shape  of  breaking  crest  length  distributions  of  the  two 
data  sets,  likely  associated  with  different  development  stages  of  the  wave  field.  We  calculated  overall 
mean  distributions  from  24  SBC  runs  with  Cp/u*  <  25,  characterizing  growing  seas,  and  from  14  HI 
runs  with  cp/u*  >  45,  characterizing  mature  seas.  These  are  shown  in  the  left  panel  of  Fig. 3.  There 
were  no  comparably  mature  sea  states  from  the  SBC  data  set  to  include.  In  growing  seas,  breaking 
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spans  the  entire  spectrum  from  dominant  waves  to  scales  associated  with  0.  lcp,  where  cp  is  the  speed 
of  the  spectral  peak  waves,  whereas  in  mature  seas  no  breaking  is  observed  at  scales  larger  than  those 
corresponding  to  0.6cp. 

The  normalized  fifth  moment  of  the  breaking  crest  length  distributions  is  shown  in  the  right  panel  of 
Fig.  3.  The  fifth  moment  of  the  A(Cb)  distribution  is  related  to  the  spectral  energy  dissipation  rate.  The 
strongest  dissipation  rate  is  seen  to  occur  at  small-to-intermediate  scales  in  the  mature  sea,  but  at  the 
large  wave  scale  in  the  developing  sea.  Thus  the  velocity  scale  of  the  peak  of  the  breaking  dissipation 
rate  decreases  with  increasing  wave  age.  Despite  the  breaking  in  mature  seas  occurring  over  a  smaller 
spectral  bandwidth  than  in  growing  seas,  the  total  dissipation  rate  in  the  mature  seas  is  about  three 
times  higher.  This  is  attributed  to  higher  energy  and  momentum  fluxes  from  the  wind  in  HI. 
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Figure  3.  Overall  mean  breaking  crest  length  distributions  (left panel)  and  their  5th  moment  (right 
panel),  for  subsets  of  the  data  representing  developing  seas  (gray  triangles)  and  mature  seas 

(black  circles). 
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Figure  4.  Bin-averaged  exponent  m  (slope)  of  the  breaking  crest  length  distribution  as  function  of 
wave  age  ( c/u *)  for  the  Santa  Barbara  channel  (gray  triangles)  and  Hawaii  experiment  (black 
circles).  The  dashed  line  depicts  the  classical  Phillips  [1985]  result  m=-6.  The  vertical  bars  indicate 

± 1  standard  deviation. 


Further,  in  Fig.  4,  it  is  seen  that  the  breaking  crest  length  spectrum  tends  to  fall  off  more  slowly  than 
the  c'6  behavior  predicted  by  Phillips  [1985].  This  result  was  statistically  significant  at  the  95% 
confidence  level.  We  note  that  this  conclusion  is  weakly  sensitive  to  the  data  processing  methodology 
used,  as  discussed  below. 

Breaking  strength  parameter  b 

We  note  that  the  transition  from  the  kinematic  quantity  of  the  crest  length  distribution  to  the  dynamics 
of  energy  dissipation  rate  requires  knowledge  of  the  breaking  strength  parameter  b.  Estimates  of  the 
spectrally-resolved  breaking  parameter  are  only  starting  to  become  available  and  all  previous  studies 
based  on  field  data  report  the  scale-integrated  value  beff.  Reported  values  of  befffrom  different 
experiments  span  more  than  3  orders  of  magnitude  between  these  experiments,  from  3x10'  to  7x10'“ 

[ Gemmrich  et  al.,  2008;  Melville  and Matusov,  2002;  Thomson  et  al.,  2009].  Recent  work  by  Romero 
et  al.  [2012]  suggests  that  approximating  b(Cb)  by  a  unifonn  level,  independent  of  Cb,  is  a  reasonable 
approach,  at  least  for  wave  scales  shorter  than  the  dominant  waves.  In  any  event,  from  our  direct 
measurements  of  A(cb)  and  surface-piercing  dissipation  rates  in  the  mature  wind  seas  off  HI,  we 
directly  measured  the  effective  breaking  strength  parameter  beff  to  be  4.2  ±1.8  x  I  O'"  for  the  wave  age 
(Cp/u*)  range  of  40-47. 

Sensitivity  to  data  analysis  methodology 

In  view  of  the  rather  large  scatter  in  beff  reported  by  different  investigators  pointed  out  above,  we 
revisited  the  sensitivity  of  the  whitecap  analysis  methodology.  While  there  are  differences  in  the 
detailed  measurement  of  the  individual  breaker  geometry,  a  significant  source  of  variation  arises 
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because  some  investigators  have  chosen  to  partition  the  time-dependent  geometry  of  each  breaker  into 
variable  speed  bins  as  the  breaker  slows  down,  while  other  investigators,  including  ourselves,  assign 
the  mean  crest  length  of  a  given  breaker  to  its  initial  speed,  as  envisaged  by  Phillips  [1985].  The 
variable  speed  allocation  method  has  the  effect  of  aliasing  A(Cb)  to  slower  waves,  thereby  increasing 
the  exponent  -m.  Our  analysis  based  on  the  RaDyO  data  sets  estimates  this  difference  at  0(0.3  to  0.4). 
The  corresponding  impact  on  befr  for  the  RaDyO  data  was  a  factor  of  0(2)  increase  associated  with  the 
variable  speed  method.  We  are  investigating  other  possible  sources  that  may  be  contributing  to  the 
large  dynamic  range  of  reported  beff  values. 

IMPACT/APPLICATIONS 

This  effort  will  contribute  a  more  detailed  characterization  of  the  wind  driven  air-sea  interface, 
including  wave  breaking  (whitecaps  and  microscale  breaking).  This  is  needed  to  provide  more 
complete  parameterizations  of  these  processes,  which  will  improve  the  accuracy  of  ocean  optical 
radiative  transfer  models  and  trans-interfacial  image  reconstruction  techniques. 
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